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Abstract: We analyze the 1+1 CFT states dual to hot (time-symmetric) 2+1 multibound¬ 
ary AdS wormholes. These are black hole geometries with high local temperature, u > 1 
asymptotically-AdSs regions, and arbitrary internal topology. The dual state at t = 0 is dehned 
on n circles. We show these to be well-described by sewing together tensor networks corresponding 
to thermoheld double states. As a result, the entanglement is spatially localized and bipartite: 
away from particular boundary points (“vertices”) any small connected region A of the boundary 
CFT is entangled only with another small connected region B, where B may he on a different 
circle or may be a different part of the same circle. We focus on the pair-of-pants case, from 
which more general cases may be constructed. We also discuss hnite-temperature corrections, 
where we note that the states involve a code subspace in each circle. 
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1 Introduction 

The thermofield double (TFD) state 

\TFD) = (1.1) 

E 

on two copies of a quantum field theory serves as the poster child for many ideas [1-7] relating the 
emergence of bulk geometry to entanglement in some dual theory. As explained in [1] , although a 
single copy of a CFT can be naturally dual to bulk quantum gravity with a single asymptotically 
AdS boundary, the particular entanglements between the two copies described by (1.1) allow 
it to be dual to a two-sided eternal black hole in which two distinct asymptotic regions are 
connected by an Einstein-Rosen bridge^. The state also typifies relations between the area of 
codimension-2 surfaces and CFT entanglement encapsulated in the Ryu-Takayangi conjecture [7] 
and the covariant generalization by Hubeny, Rangamani, and Takayanagi (HRT) [5]. Here and 

^Though there may be interesting subtleties; see e.g. [8]. 
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below we work in the regime where the bulk planck scale Ip is small in comparison with the bulk 
AdS scale EacIS (which we generally set to 1), or equivalently where iV 3> 1 in the CFT (i.e., large 
central charge c for a 1 + 1 CFT). 

In discussing \TFD), it is natural to focus on the bipartite entanglement between the 
associated two copies of the CFT. This entanglement has a special structure: as shown in [9], the 
entanglement is both local and bipartite in the sense that, when studying regions of the CFT 
of size greater than the thermal scale, a given region can be said to be entangled only with the 
corresponding spatial region in the second CFT. In particular, when we consider regions A, B (in 
the same or opposite CFTs) separated by more than this scale, the mutual information 

J(A : B) = S{A) + S{B) - S{AB), (1.2) 

vanishes at leading order in large N. This result can easily be understood from a CFT path 
integral point of view. In general, the thermofield double state is calculated by a CFT path 
integral over a cylinder, linking the two copies of the spatial section the state (1.1) is defined 
on. In the high temperature limit, this cylinder becomes short compared to its circumference, so 
when we consider regions larger than the length of the cylinder, the resulting state naturally only 
entangles regions on one boundary with the corresponding region on the other boundary. 

It will be useful below to visualize this result in the language of tensor networks; see e.g. 
[10]. The rather trivial nature of the above entanglement then translates into a similarly-trivial 
coarse-grained tensor network description of \TFD) as shown in figure 1. 

While the thermofield double state is a useful simple example, it is important to find further 
examples where we can understand the relation of bulk geometry to CFT entanglement structures. 
We are also interested in exploring the role played by multi-party entanglement in connections 
between 3 or more subsystems and what form it takes in the associated CFT states, see e.g. 

The vast literature on holographic entanglement has focused primarily on bipartite relations 
between a given subsystem in the CFT and its complement, so that relatively little is known about 
multiparty issues. One general result is the monogamy of holographic entanglement established 
in [14]. But a more detailed investigation of multipartite entanglement was recently initiated 
in [12] using a class of 2-|-l-dimensional black hole spacetimes [15-18] describing a collapsing 
wormhole that connects n regions each asymptotic to (global) AdSs. When the corresponding 
Euclidean geometries define the dominant saddle of a natural path integral, such geometries 
are dual to entangled states on n copies of a 1-1-1 dimensional CFT on x M dehned by a 
path integral on a Riemann surface S with n circular boundaries [1, 19-21]. The corresponding 
entanglement was found to display a rich dependence on the moduli, including regimes of purely 
bipartite entanglement, and others of strong multipartite entanglement. Interestingly, the strongly 
multipartite regions identified in [12] corresponded to bulk black holes with temperature less than 
the AdS scale^. The recent work [22] describes an infinite family of generalizations of results from 

^It remains an open question whether such phases ever dominate the path integrals described above. But even if 
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Figure 1: A simple tensor network displaying the localized pnrely-bipartite entanglement 
characteristic of holographic \TFD) states at large N on scales longer than the thermal scale. 
Each node represents a region in the CFT of scale longer than the thermal scale. We focus mainly 
on CFT states on x M where one takes a high-temperature limit in order to ht many such 
long-distance regions onto the circle, though one may equally-well consider the planar case. The 
solid links are the entangling tensors implied by (1.1). The dashed lines guide the eye by linking 
neighbouring regions in each of the two CFTs. 


both [14] and [12]. 

We focus below on the opposite limit in which all bulk black holes have high temperature. 
The length of their horizons is then very large with respect to the AdS scale. We will show that 
these geometries are dual to states constructed by sewing together copies of \TFD), as shown in 
figure 2. The entanglement is thus both local and bipartite away from small regions containing 
certain “vertices” where the sewing involves three or more copies of \TFD). From the CFT path 
integral point of view, this arises because the boundary circles are large compared to the distance 
between them; in a conformal frame where the boundaries are hnite size, there are thin strips 
joining them, corresponding to the short tensor networks in figure 2. In section 4 we will justify 
this picture more quantitatively by showing that local pieces of the surface S are described by 
regions of BTZ up to exponential corrections. As a result, as in [23] tripartite entanglement 
appears to localise in isolated AdS-scale regions of the bulk. Away from these vertices, the 
construction of the state involves only the sewing together \TFD)''s of inverse temperature /3i 
and (32, giving a local version of the \TFD) of inverse temperature (3i + (32- Since we focus on 
l-fl CFTs, we henceforth refer to the limit of large central charge c rather than large N. 

Note that nothing prevents sewing operations that link together disjoint regions in the same 
CFT as shown in figure 2 (bottom). As we will see, this also provides an interesting picture in 
our limit of CFT states dual to single-boundary black holes with internal topology. The reader 
should thus be aware that, while we use term “multiboundary” below, this explicitly includes the 
very interesting case n = 1 as well as n > 2. 

One may expect each local piece of \TFD) in figure 2 (right) to correspond to a bulk region 

not, one presumes them to be dual to some other class of CFT states whose entanglement must be correspondingly 
multipartite. 
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Figure 2: Two topologically-distinct ways in which three copies of the tensor network in figure 1 
can be sewn together (left figures) into a single tensor network (right figures) defining a state 
on 3 copies of the system. The dashed lines (red in color version) internal to the left diagrams 
guide the eye toward recognising the 3 constituent copies of the network in figure 1. Links that 
meet across adjoining pairs of dashed lines are contracted, establishing entanglement between the 
remaining boundaries (marked 1, 2, and 3). In the bottom-left figure, two parts of the outermost 
tensor network are contracted with each other, resulting in two well-separated regions of boundary 
1 becoming entangled with each other as shown in the bottom-right figure. As discussed below, 
all 3-boundary time-symmetric vacuum wormholes with pair-of-pants topology (orientable with 
no handles) and large horizons correspond at the moment of time-symmetry to one of the cases 
shown, or to the degenerate case that interpolates between them, when described in the “round” 
conformal frame in which the energy density is taken to be constant along each of the 3 boundaries. 
Although we show only a simplified cartoon of the full tensor network, we argue below that sewing 
the actual \TFD) tensor networks together in this way describes the corresponding CFT states 
with exponential accuracy away from the two ‘vertices’ in each diagram where 3 \TFDys meet. 


whose geometry near t = 0 is well-approximated by a corresponding piece of BTZ. We show in 
section 5 that this is indeed the case, and thus that bulk Ryu-Takayanagi (or, more precisely, 
HRT) calculations are consistent with the entanglements shown. 

We begin by reviewing aspects of general multiboundary wormholes and their relation to 
CFT states in section 2. Section 3 then studies the high temperature (equivalently, large horizon 
length L) limit of the geometry of S. We show that the region between the horizons in S becomes 
unimportant in this limit. This allows us to argue in section 4 that the CFT path integral 
produces the structure described by figure 2. Section 5 then describes how this same result is seen 
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in the bulk HRT calculation and argues that the desired bulk wormhole does indeed dominate the 
corresponding bulk path integral. Section 6 briefly addresses flnite temperature corrections and 
we conclude in section 7 with discussions of the general n-boundary case, internal topology, higher 
dimensions, and future directions. In particular we comment explicitly on examples with n = 1. 


2 Path integrals, states, and bulk geometries 

We now commence our review. As is well known, the thermofleld double state of inverse 
temperature (3 is computed by the CFT path integral on the cylinder of circumference 27r and 
height (3/2. Here one regards each of the two circular boundaries as the (say) t = 0 slice of a 
corresponding CFT. The path integral between held configurations 0i,02 on the two boundaries 
then gives the wavefunction ^(^i, 02) of the joint state of the two CFTs.^ At sufficiently high 
temperatures, the corresponding bulk path integral is dominated by a saddle point associated 
with the Euclidean BTZ black hole. In this case we may say that, to good approximation, the 
corresponding Lorentz-signature bulk black hole is dual to \TFD). 

The cylinder of circumference 27r and height 0/2 plays two important roles in the BTZ 
geometry. First, it is conformally equivalent to (half of) the boundary of Euclidean BTZ. This 
is what allows Euclidean BTZ to be a saddle for the desired bulk path integral. But this same 
cylinder is also conformal to the BTZ geometry at t = 0, which may be equally-well considered 
as a slice of either the Euclidean or the Lorentzian black hole. This may be seen by recalling 
[24] that Euclidean BTZ can be constructed as a quotient of global Euclidean AdSa (i.e., of the 
hyperbolic three-space by an isometry. The simplest statement requires two steps. One first 
writes Euclidean AdSs in terms of its slicing by hyperbolic planes (equivalently, by copies of 
Euclidean AdS 2 ) as 

^2 — = -I-cosh^ (2.1) 

^AdS 

where is the metric on the unit-radius H^. One then quotients each slice by a discrete 
group F = : n G Z} generated by some hyperbolic element g of its S'L(2,M) group. ^ The 

action of g and its fundamental domain in are indicated in figure 3. Since the different 
slices in (2.1) differ only by the overall scale factor cosh^t®, the same is true of their quotients. 
The spatial slice at t = 0 (equivalently, Ie = 0) is thus conformal to the geometry at = —oo. 
This is half of the Euclidean boundary, with the other half being Ie = -l-cxo. We may therefore 
write \TFD) as given by the CFT path integral over the Riemann surface defined by the t = 0 
slice of the corresponding BTZ geometry. 

^In this discussion we assume for simplicity that the CFTs admit (anti-unitary) time-reversal symmetries T 
which can be used to map bra-vectors to ket-vectors and vice versa, and which can therefore be used to construct 
(1.1) from the thermal operator 

“^In other words, thinking of SL{2,M.) as the Lorentz group SO{2, 1) of 2+1 Minkowski space this g must be a 
boost preserving some spacelike direction. 
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Figure 3: The quotient of the hyperbolic plane by T. The pair of labeled geodesics are 
identified by so the region between them forms a fnndamental domain for the qnotient. The 
minimal closed geodesic H is the horizon for the resnlting BTZ geometry. 

This hnal conclnsion can be extended to a mnch larger class of states. Any Riemann snrface 
S with n boundaries can be written as a quotient of by some discrete subgroup of SL{2, M). 
We may use (2.1) to construct a corresponding quotient of Euclidean AdSs, with S conformal 
to both the slices at f = 0 and Ie = —oo. So long as this saddle point dominates the bulk 
path integral with boundary conditions dehned by the = — oo slice, to good approximation 
the corresponding Lorentz-signature bulk solution - given by substituting = —it into (2.1) 
" is dual to the CFT state dehned by the path integral over the slice at f = 0. For notational 
simplicity we identify S with this slice below and write the CFT state as |S). These qnotients of 
AdSa were hrst considered in [16], and the holographic relation to |S) was introdnced in [19-21]. 
An exploration of the entanglement properties of these states was initiated in [12]. 

The Lorentz-signature solutions describe wormholes connecting n asymptotically-AdSa bound¬ 
aries. By topological censorship [25, 26], each boundary is associated with a distinct event horizon. 
A special property of AdSa vacuum solutions is that the geometry outside each event horizon 
is precisely that external to some BTZ black hole. This allows us to dehne a useful “round” 
conformal frame, in which the usual rotational symmetry of this BTZ region is a symmetry of 
the boundary. That is, for each of these exterior regions there is a coordinate 0* such that 
is an exact rotational Killing held in the region outside the horizon (and in fact in an open 
neighbonrhood in the interior of the horizon as well). The ronnd conformal frame is the one in 
which the CFT lives on a spacetime with standard cylinder metric 

ds^ = —dt^ -I- d(j)f (2.2) 

with (j) ^ (j) + 271. In addition, the BTZ exterior implies that the bifurcation surface of each 
horizon - where the fntnre and past horizons meet - is a geodesic in the t = 0 snrface. The key 
novelty in the n > 2 cases is the existence of a “cansal shadow” region in between these horizons. 
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Figure 4: The surface S as a quotient of the Poincare disc for n = 3. The pairs of labeled 
geodesics (blue and red in colour version) are identihed by the action of T. The region of the 
Poincare disc bounded by these geodesics provides a fnndamental domain for the qnotient. B^, 
B 2 and Bi = i?i+ U i?i_ become the desired three circnlar bonndaries. There are corresponding 
minimal closed geodesics H 3 , H 2 and Hi = U ifi_. The lengths La of these geodesics fully 
characterize the geometry of S. 

Our ideas will apply to a codimension one limit in the moduli space of such Riemann surfaces 
for any n, but for simplicity we will focus our discussion on the case where E is an orientable 
surface with three boundaries and no handles. Such surfaces are topologically the same as a 
pair of pants. This is the simplest non-trivial example, and is also a primitive building block 
for constructing other cases, since a general orientable Riemann surface can be constructed by 
sewing together pairs of pants. The relevant qnotient of is depicted in hgnre 4. The modnli 
space of pair-of-pants Riemann surfaces can be parametrized by the lengths (a = 1, 2, 3) of 
the three horizons, which as nsnal we take to be measnred in nnits with = 1- Withont loss 
of generality we take L 3 > Li, L 2 . The cansal shadow is the region in between these geodesics. 

Properties of such states were explored in [12], with most emphasis on the so-called puncture 
limit La ^ 1. In particular, [12] showed that in this limit E is conformal to the Riemann sphere 
with small holes removed around n points, and hence |E) can be related to an n-point function in 
the CFT. For the three-bonndary case, the state was determined np to some constant factors to 
be, in the ronnd conformal frame, 

|S) = (2.3) 

ijk 

where 

/3a = /3a-21nrrf-21n3, (2.4) 

and Cijk are the three-point OPE coefficients, /Sa = ^ is the inverse temperature of the BTZ 
geometry associated with the region near the a’th bonndary, and is an nndetermined constant 
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independent of the moduli. The rather explicit expression (2.3) exhibits both dependence on the 
structure of the CFT and Boltzmann-like suppression factors similar to the thermoheld double 
state. 

3 Geometry of E in the high temperature limit 

Our current aim is to elucidate the structure of |S) in the limit —)■ cx) with fixed ratios La/L^. 
This is the opposite of the limit emphasized in [12, 23]. We assume L 3 > Li,L 2 , so the ratios 
L 1 /L 3 , T 2 /T 3 take values in (0,1]. In this limit, the geometry of S again simplihes. The essential 
point is that the causal shadow region will play a relatively unimportant role. We will focus on 
the pair of pants case, but the discussion is easily extended to arbitrary Riemann surfaces. We 
comment on this extension in section 7. 

Our limit can be characterised as a high temperature limit, in the sense that the BTZ horizon 
in each of the exterior regions becomes large compared to the AdS scale (as for a high T BTZ 
black hole). But we note that the restriction of the state |S) to a single boundary is not necessarily 
even approximately thermal: as discussed in [ 12 ], when one La is larger than the sum of the other 
ones, reduced density matrix in that exterior region has much less entropy than the thermal value 
at the same energy. 

To understand the geometry of S in our limit, it is useful to introduce a different presentation 
using two patches with BTZ coordinates on each^. We split hgure 4 in half along the horizontal 
geodesic (not drawn explicitly) joining boundaries Bi and B 2 . This divides S into two identical 
regions E±, each containing half of each horizon Ha. The surface E is then recovered by gluing 
together E± along three geodesics, the two identified geodesics in figure 4 and the new split. We 
label these geodesics Gab = Gba with a 7 ^ 6 labelling the boundaries they run between; see hgure 
5 (left). They can be described more formally as the hxed points of a Z 2 isometry of E, which 
acts as a rehection 0 —)■ 27r — 0 in the round conformal frame on each of the boundaries (with 
appropriate choices of the origin 0 = 0 on each boundary). The event horizon Ha of boundary a 
is the unique geodesic that runs orthogonally between the two geodesics Gab, Gac ip p c) that 
end on boundary a. Our partition of E into E± also breaks each horizon Ha into two pieces Ha±. 

It is useful to describe E± in planar BTZ coordinates.® Consider for dehniteness E+. We 
choose the BTZ coordinates to be 

^AdS w i 

with p G (—C) 0 , 00 ). Thus our reference BTZ solution has inverse temperature 27r. Without loss of 
generality, we take L 3 > L 2 > Li and orient E_|_ such that the portion of H^ in E+ lies along the 

®In the actual history of our project, this description was also inspired by computing mutual information on 
pairs of pants with large La using the technology of [27]. 

®By planar BTZ we mean the usual BTZ coordinates with no identification on the spatial coordinate on the 
boundary; this provides a coordinate system on the whole of thought of as the t = 0 surface in AdSa. Since 
our Y,± are subregions of they can be conveniently described in these coordinates. 






Figure 5: The region S+ bounded by the geodesics Gab, half each of B 2 , - 83 , and i?i+ shown in 
the Poincare disc (left) and the BTZ frame strip (right). The BTZ presentation is chosen to place 
the half-horizon along the BTZ horizon. The geodesics G 13 , G 23 are respectively the lines 
X = X = ^. In contrast, G 12 lies in the upper half of the strip; its endpoints have x = Xi,X2 
with p = -|-oo. Half each of Bi, B 2 is mapped respectively to the line segments x G 
X G [x 2 , at p = cx) respectively, whilst half of B 3 is mapped to p = — 00 . The corresponding E_ 
is the symmetric region below G 23 in the Poinare disk (left) and has an identical representation 
in the BTZ strip. 


horizon at p = 0, and the boundary B 3 lies at p = — 00 , in both cases for x G [— L3/4, T3/4]. Since 
the geodesics G 13 and G 23 intersect H 3 orthogonally, they will lie at x = —L 3 /A and x = L 3/4 in 
these coordinates. The other two boundaries Bi and B 2 lie at p = cxd, for x G [—L3/4, xi] and 
X G [x 2 ,T 3 / 4 ] (with xi < X2), and the remaining geodesic G12 runs between these points xi,X2. 
The portions of Hi, H 2 in are the geodesics running from the edges of the strip to meet G 12 
orthogonally. These coordinates are illustrated in hgure 5. 

The half-surface E+ is thus a strip x G [—L3/4, L3/4] in the planar BTZ coordinates, with 
a bite cut out of the middle above G 12 . It is important to emphasize that the boundaries at 
—L3/4, L3/4 are not identihed with each other; instead they and G 12 are identihed with the 
corresponding boundaries in a second copy of this region. 

As we verify in appendix A, varying the endpoints Xi,X 2 of i?i,i ?2 generates all possible 
lengths Li, L 2 for the remaining horizons Hi, H 2 and the map (xi, X 2 ) e->- (Li, L 2 ) is both smooth 
and one-to-one. When we take the limit of large La (at hxed ratios), the results simplify, with a 
form that depends on the relative lengths. For L 3 — {Li + L 2 ) 3> 1, 

~ y - Y X 2 ~ Y “ y+ ^°S2, (3.2) 

where the tildes (~) represent agreement up to exponentially small corrections. Note that up to 
a fixed order-one offset, the endpoints are respectively Li/2 and L 2/2 from the ends of the strip. 
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(a) Li + L2 > 

Xi X 2 



(b) Li + L2 ~ -bs 



(c) L3 > Li + L2 

Figure 6: Half of the pair-of-pants (either S+ or S_) described as a region in planar BTZ. 
Three examples are shown representing distinct regions of modnli space: Li + L 2 > (top), 
Li + L 2 ~ T 3 (middle), L 3 > Li + L 2 (bottom). 

In the complementary case L 3 — (Li + L 2 ) <C 1, we find instead 

La T L2 — T3 \ , . 

^- f - -j . (3-3) 

In our BTZ presentation, the long length of corresponds directly to the large width of the 
strip. The horizons Hi, H 2 are also long as a result of extending over a large coordinate distance 
in the x direction and possibly also from extending out towards the boundary of the strip at 
large p. If both of them together are shorter than (Li + L 2 < T 3 ), they terminate on Gu 
in the interior of the strip, staying within an order one distance from the horizon at p = 0 
along their whole length, as in the last panel of hgure 6 . When the sum is larger (which includes 
the case where the three lengths are equal), the length of the interval X 2 — Xi is exponentially 
short, hence Hi, H 2 meet G 12 at large p, as in the hrst panel of hgure 6. In both cases. Hi and 
H 2 approach H^ exponentially for \x — xi^ 2 \ ^ 1 - 

The contributions to Li, L 2 from the width of the strip or from Hi and H 2 running to large 
p look different, but we should remember that treats the three horizons symmetrically, so 


Xi +X2 Li - L2 X2 - Xi 


= exp 
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this is just an artifact of our choice of coordinates. The symmetry can be made manifest in an 
appropriate Poincare disk representation; see e.g. hgnre 8 of [12], The relationship between any 
pair of horizons is thus much the same; consider for example Hi and H^. We can see explicitly 
from the calculation in appendix A that the minimal distance between them is exponentially 
small, and that they remain exponentially close over a large region. Thns, the area of the cansal 
shadow region remains hnite even as their length becomes large. 

In fact, since the boundaries of the causal shadow are closed geodesics (and thus have vanishing 
extrinsic cnrvatnre), the Ganss-Bonnet theorem reqnires this area Acs to be independent of the 
moduli La (for any hxed genus g and number of boundaries n). For the pair of pants we hnd 
Acs = 27r; more generally Acs = 2 (n — 2 + 2 g) 7 T. As we will see in the next section, this implies 
that the cansal shadow region plays little role in the path integral constrnction of the CFT state 
|S). 

In the case where L 3 — (Li + L 2 ) 3> 1, the endpoints of the geodesic G 12 are far apart in 
coordinate distance, and it will also be exponentially close to over most of its length. When we 
glue S+ and S_ to form S, the section of H 3 that is close to G 12 will lie close to the corresponding 
section of in S_, as in hgure 2 (right). All remaining cases with L 3 > L 2 > Li are intermediate 
between the two jnst described. 

4 The CFT state at large La 

Let us now consider the implications of the above results on the structure of S for the CFT state 
|S). In this section we will argne for large La that |S) will be described to exponential accnracy 
by hgnre 2 (right). In particnlar, when restricted to appropriate regions it agrees to exponential 
accnracy with the corresponding restriction of a thermoheld donble state \TFD). We also show, 
under the assumption that non-handlebody contributions can be ignored, that the Euclidean bulk 
geometry (2.1) dominates the bnlk path integral dehned by nsing E as the conformal bonndary. 
It follows that, to exponential accuracy, our bulk pair-of-pants wormhole is dual to the state 
described by hgure 2 (right). 

Recall that |S) is dehned by the CFT path integral over S. We will use the BTZ representation 
of E± associated with hgnre 6 to break E into three pieces Ei^ 2,3 that are topologically cylinders, 
corresponding to hgnre 2 (left). Each piece E^ will contain the entire pair-of-pants bonndary Ba, 
bnt no portion of the bonndaries Bf, for b ^ a. The decomposition is dehned by drawing a graph 
in E as shown in hgure 7. As noted in the hgure caption, far from the vertices the piece La nearly 
coincides with the region Ea of E exterior to Ha (i.e., lying between Ha and boundary a). The 
geometry in this latter region is just that of the appropriate BTZ solution outside the horizon 
and is conformal to a ronnd (rotationally-invariant) cylinder. 

We wish to regard both Ea and E^ as path integrals constructing states \Ea), |E(j), each 
of which is dehned on two copies of our CFT (on the outer and inner boundaries of Ea or E^ 
respectively). Indeed, we may write |Ea) = Sa\Ea) where Sa is the operator dehned by the path 
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Figure 7; The decomposition of into pieces Ei±, Ei±, Ei±. We then glue Ea+ to Ea_ along 
the relevant Gab to make pieces E^ conformal to cylinders with E^ containing all of boundary a. 
The decomposition is determined by a graph. In cases (a) and (c), the graph has two trivalent 
vertices. In case (a) each piece E± contains one vertex v±. The 3 edges of the graph each connect 
to V- running between two distinct horizons Ha,Hb for a ^ b. In case (c) both vertices v,v' 
lie on the cut along Gu, as does the edge that connects them. The other two edges are loops 
connecting u to u or v' to v'. One lies between Hi and while the other lies between H 2 and 
H^. Case (b) represents a degenerate limit interpolating between the two in which we choose to 
fuse the two vertices into a single 4-valent vertex lying on G 12 . The graph has two edges, each 
of which are loops. One lies between Hi and H^ while the other lies between H 2 and H^. In 
the regions far from the vertices the pieces Ei 2,3 dehned by cutting E along the edges of the 
appropriate graph differ from the cylinders dehned by the regions outside horizons i 7 i, 2,3 only by 
exponentially small amounts. 


integral over the causal shadow region Sa = T,a/Ea in E^ beyond the horizon Ha- We specify 
the conformal frames of all states by again taking (pa to dehne the standard angle on the CFT 
cylinder; this involves a natural extension of (pa through the causal shadow Sa- The region Sa 
is topologically an annulus and so can be conformally transformed to a cylinder. But Sa is 
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exponentially thin over most of its circumference; indeed, setting a = 3 (so that we may replace 
03 by 27r/L3 times the BTZ x) and multiplying the BTZ metric (3.1) by ^aL( 1 +p") Ogives a 
metric 

dsls = (4.1) 

where y = tan“^ p ranges over [ 0 , f{x)] with f{x) is exponentially small far from the vertices of 
our graph. Introducing y = y/f along with x such that dx = dx/f, and multiplying (4.1) by 
gives a metric 

r‘^ds% = dx^ + {dy + f'ydxf, (4.2) 

where /' = df /dx, on cylinder of unit height y G [ 0 , 1 ] but with exponentially large circumference. 
The metric is not flat, though it differs from the standard cartesian flat metric dx^ + d/p only by 
exponentially small corrections proportional to powers of /'. It follows that there is a further 
conformal transformation to a metric cylinder of unit height - and with exponentially small 
difference in circumference from the range of x - whose action on the region far from the vertices 
is exponentially close to the identity map^. 

Rescaling this cylinder to one of circumference 27r allows us to write the path integral over Sa 
in terms of the operator with exponentially small /3. Up to exponentially small corrections, 
this operator acts as the identity with respect to degrees of freedom associated with spatial regions 
of order-one size as measured by the original spatial coordinate <pa- So far from the vertices we 
may identify \Ea) and |Sa) with exponential accuracy. 

On the other hand, the exterior region Ea is half of the BTZ t = 0 surface, so \Ea) is given 
by a path integral over a cylinder of length 0/4, so it is a copy \TED)a of ( 1 . 1 ) at twice the 
temperature of the associated bulk horizon. Thus |Sa) is, up to exponentially small corrections, a 
thermoheld double state. Recall that for large La such \TED)a are described by a tensor network 
of the form shown in £g. 1 . 

It remains only to sew the |Sa) = \TED)a together into |S). The sewing procedure is dehned 
by the way the path integrals Si 2,3 combine to form S. Away from the vertices of the graph, this 
identihes the horizons in £' 1 ^ 2 , 3 : for Li -|- L 2 > R 3 , parts of H 3 are identihed with each of Hi and 
£ 2 , and the remaining parts of Hi and H 2 are identihed with each other. For Li -|- L 2 < R 3 , Hi 
and H 2 are each entirely identihed with corresponding parts of £ 3 , and the remaining regions 
of £3 along Gi 2 are identihed with each other. Since the sewing operation on path integrals 
coincides with the sewing operation on tensor networks - one simply sets all arguments equal 
along the seam and integrates over allowed values®- this implies that the state |E) is given to 
exponential accuracy by hg. 2 (right), with the top picture relevant for Li -|- L 2 > R 3 and the 
bottom picture relevant for Li + L 2 < L^. 

^Here we use the fact that conformal transformations satisfy an elliptic equation with a Green’s function that 
decays exponentially along a strip. We expect that similar arguments are common in the literature, but for a 
specific example the interested reader interested in details may consult for comparison e.g. section 3.1.1 of [12]. 

®Here for simplicity we again make use of the time-reversal symmetry mentioned in section 2 to turn bra-vectors 
into ket-vectors. 
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Finally, we also wanted to see that |S) is dnal to onr bnlk geometry with moment of time 
symmetry E. As in [12] we assnme that the dominant saddle of the associated bulk integral 
is a handlebody. The other possible bulk saddles discussed in [12] correspond to disconnected 
Lorentzian geometries. It is natural to expect this saddle to dominate at large temperatures, 
by analogy to the familiar result for \TFD) that disconnected solutions dominate only at low 
temperatures. But one can now make a further argument based on entanglement. If the HRT 
proposal is correct, and in particular if entanglement is associated with extremal surfaces in the 
real Lorentz-signature geometry, the disconnected geometries cannot reproduce the entanglement 
structure of hgure 2 (right), which involves entanglement between the different boundaries at 
leading order in the central charge. It would be interesting to verify this conclusion by direct 
computation of the Euclidean actions, as it would serve as a check on HRT. 


5 Holographic entanglement calculations 

The previous section used the GET path integral to show that the GET state |S) is given by 
figure 2 (right), so that the state has local bipartite entanglement with the same local structure 
as the thermofield double state. In this section we will buttress that argument by showing that 
our picture of the geometry of E, now thought of as the t = 0 surface in the bulk spacetime, 
gives consistent results for entanglement from holographic Ryu-Takayanagi calculations. Indeed, 
in the history of our project we originally discovered that the state had this simple bipartite 
structure by performing these holographic calculations explicitly. We consider the entanglement 
for a region in boundary 3, since our coordinates are adapted to this boundary, but by symmetry 
similar results apply in the other cases. 

Gonsider first a region in boundary 3 where the horizon is exponentially close to either Hi 
or H 2 , that is Xi — a; 1 or a; — a ;2 f. In the exterior region E^, the planar BTZ coordinate 
X is identical (up to a scale and a shift of origin) to the round conformal frame coordinate 03 
dehned in section 2: 03 = |^a:. In the other exterior region, at similar x but outside Hi, H 2 , 
because the horizon H 12 is exponentially close to H^, the planar BTZ coordinate x agrees with 
012 (up to a scale and a shift of origin) up to exponentially small corrections. We can take for 
example 0i = j^x. This is manifestly true near p = 0 (see (3.1)) and continues to hold at large p 
due to the properties of geodesics in hyperbolic geometry®. So for xi — a; S> 1, x — X 2 3> 1 we 
may take x to define the round conformal frame on all three boundaries up to exponentially small 
corrections. 

®Two geodesics on fired at slightly different angles from the same point will diverge exponentially as 
measured by the proper distance separating them as the curves approach the boundary. So curves of constant 
01^2 and X that meet at the horizon also diverge in a similar manner near the boundary. But two geodesics hred 
orthogonally from different points x', x" of the horizon again diverge exponentially at precisely the same rate. So 
a curve of constant 0 i _2 that meets the horizon at x' with \x' — xi ^2 \ ^ 1 will meet the boundary at some x” with 
\x'' — x'\ still exponentially small. 


14 



The above relations allow us to easily map those geodesics involved in any HRT calculation 
of the mutual information between subregions of boundaries 1 and 3 (or 2 and 3) that lie far from 
Xi,X 2 to geodesics in BTZ. The BTZ calculation was studied in [9], who found that for regions 
much larger than the thermal scale, the mutual information is simply proportional to the size of 
the overlap between the two regions. The overlap is maximal when the two regions are directly 
opposite each other, in which case the high-temperature result (3.27) of [9] becomes 

I(A : B) = S(A) + S(B) - S(/l U B) = ^ - A.^)) ^ ^ 

Applying appropriate scalings to (5.1), the mutual information between corresponding regions of 
boundaries 1 and 3 with xi — a; S> 1 is 

HA : B) = ^Ax + 0(1) = + 0(1) = + 0(1). (5.2) 

In addition, since the region of boundary 1 with xi — x S> 1 is well-separated in the bulk from 
the region of boundary 2 with x — X 2 3> 1 , it also follows that these two regions share no mutual 
information. The situation is exactly similar for the region in boundary 3 with x — X 2 ^ 1, which 
has a mutual information of the same form with a region in boundary 2 . 

If |x 2 — a:i| 1, there are large parts of Hi, H 2 that are far from H^, and so have yet to 
be described. This indicates that there are large intervals of 0 i ^2 along boundaries 1 and 2 
with X-values close to the endpoints xi,X 2 . But it also implies a large conformal transformation 
between the round conformal frame for Bi, B 2 and the planar BTZ coordinate x. As a result, 
the renormalized length of any geodesic connecting boundary 3 to these regions of boundaries 1, 2 
is very long and HRT calculations give no mutual information between boundary 3 and these 
regions. 

For |x 2 — Xi| 1, the above results describe the entanglement properties of boundary 3, with 
the exception of a region with length of order the thermal length scale (which in the planar BTZ 
coordinates is of order the AdS scale) near Xi,X 2 . So away from the vertices the entanglement 
structure obtained from bulk calculations corresponds precisely with that predicted by the state 
pictured in the top panel of figure 2 , given by sewing together thermofield double states. 

We now turn to the complementary case |x 2 — xi| 3> I. There is then a large region of 
boundary 3 not covered by the regions xi—x3>l,x — X 2 3>1 studied above. But across the 
region satisfying both x — Xi S> 1 and X 2 — x 3> I, the geodesic G 12 lies exponentially close to 
H 3 . So sewing together S_|_ and S_ in this region is well-approximated by simply gluing to each 
other the boundaries of (the region outside the lower half of each diagram in hgure 6 ) 
along TTsi; i.e., the result is well approximated by the region of the two-sided BTZ geometry with 

Xi < X < X 2 . 

Note that the two asymptotic boundaries of this new BTZ geometry are identihed with 
different regions of boundary 3 coming respectively from E+ and E_. In particular, since our 
slicing of the pair of pants into was performed using the Z 2 reflection symmetry, we see that 
a given value of x with xi < x < X 2 corresponds both to some point 03 and also to 27r — 03 in 
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terms of the usual coordinate on boundary 3 that dehnes the round conformal frame. Thus, in 
this case the mutual information of a region in boundary 3 with the corresponding region on the 
opposite side in boundary 3 will be as given in (5.2), supporting the local thermofield-double like 
entanglement between the two pieces of this boundary as indicated in the bottom panel of fig. 2. 


6 Finite size corrections 

We have shown that |E) has a simple structure in the high temperature limit. To use this as a 
systematic approximation to the state corresponding to hnite-size wormholes, it is interesting to 
investigate hnite-temperature corrections to this. In this section we will consider this first for the 
simple two-boundary case and then for three boundaries. 

6.1 Two boundaries 

In the two-boundary case, we want to understand and characterise the departure from the trivial 
network pictured in figure 1. The departure will be significant when we consider small regions, of 
order the thermal scale or smaller. For simplicity, we diagnose this by considering the mutual 
information between a subregion in one boundary and the whole of the other boundary. 

The key hnite temperature effect is that, for small regions, there is a competition between 
different possible minimal surfaces in the bulk homologous to A. For S(A), we need the smaller of 
/(7 a), the length of the minimal (connected) geodesic ja homotopic to A, or l(pfA'=^) + L, where 7 ^= 
is the minimal (connected) geodesic homotopic to and L is the length of the closed geodesic 
at the horizon. Similarly ^'(A'^) is determined by either /(7 Ac) or Kia) + T, and there is an 
interesting competition between these two possibilities when A is nearly the whole boundary. 

At high temperature the geodesics '^AiIa^ behave as shown in hgure 8. As a function of the 
angle 0, they drop quickly from the boundary to the horizon, hug the horizon while traversing an 
angle nearly 27r — A(/), and then quickly return to the boundary. One thus hnds 

1 (ia) = + = -f- (1 - ^) + (6.1) 

which reproduces the behaviour in (5.1) found in [9]. 

At hnite L, corrections to (6.1) are exponentially small in L when A and A'^ are larger than 
the thermal scale, and the entropy remains close to linear in A(j). But for any hnite L there is a 
Ryu-Takayanagi phase transition when either A or becomes sufficiently small. In that regime 
the relevant entropy S{A) or S'(A'^) becomes controlled by the disconnected geodesic. Thus, when 
the length of A falls below 27rlog2/L -|- 0{e~^/L) one hnds I{A : B) = 0. For A'^ smaller than 
this threshold, one hnds I{A : B) = 2S{B) = Plotting the full I{A : B) at large but hnite T 
clearly shows these “plateaux” as in hgure 9. These plateaux were studied in [28]; they can be 
characterised in terms of saturation of the Araki-Lieb inequality as discussed in [29]. 
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Figure 8: The geodesics giving, in the high-temperatnre limit, the phases of entanglement 
entropy of the nnion of a pair of intervals (cyan) lying on opposite bonndaries, along with the 
event horizon added to satisfy the homology constraint, marked by the horizontal dashed line. 
When one of the intervals is a whole boundary, there are only two relevant phases (left), otherwise 
a third phase (right) may dominate, for which the corresponding geodesics cross the horizon and 
have endpoints lying on opposite boundaries. 




Figure 9: Left: Plot for L = 10 of the mutual information /(/i; 2) between an interval Ji of size A 
in Bi and all of B 2 , as a function of A, scaled by the maximal value 2S. The mutual information 
increases approximately linearly in A and becomes non-zero at the phase transition where A ~ 
(vertical line near left edge). Right: Plot of the deviation in the mutual information shown at 
left from the high-temperature estimate This deviation is very small; for readability the 

vertical scale has been magnihed relative to the left-hand plot. The deviation is most signihcant 
for small and large values of A, and decays exponentially in L at intermediate values as expected. 


The fact that small intervals do not capture the entanglement with boundary 2 indicates 
that this information is encoded in a way that is non-local on the thermal length scale. This is 
not surprising, but it is useful to note that this correction to the large L picture has a natural 
expression in the language of [30]. It says that at hnite temperature the information about the 
entanglement between the two boundaries is not encoded locally in degrees of freedom at individual 
spatial points, but rather in a code subspace in each boundary, which entangles individual spatial 
degrees of freedom on the thermal scale. The ability to recover all of the information from any 
sufficiently large spatial subset of the degrees of freedom is the characteristic signature of such 
encoding in a code subspace. In [30], the size of the region needed to access information in a code 
subspace was related to the radial location of the bulk region encoded. Similarly, in BTZ this 
size is related to the radial position of the horizon. 

On a related note, the plateau at large A0 appears precisely when the Ryu-Takayanagi surface 
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(a) Phase 1: 1 


(b) Phase 2: 
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Figure 10: The geodesics giving the four possible phases of entanglement entropy of a single 
interval, in red, along with the event horizons added to satisfy the homology constraint, marked 
by dashed lines. 

for region A is 'j{A‘^) plus the horizon. In other words, it occurs precisely when the so-called 
entanglement wedge [3] - the region inside this Ryu-Takayanagi surface - reaches all the way to 
the horizon. Indeed, in this case we see that it touches each and every point on the horizon and 
on A’s side of the horizon it misses only a small part of the space near This suggests that the 
bulk near-horizon degrees of freedom are encoded non-locally in the CFT in such a way that they 
can be perfectly recovered from a large spatial subset A that remains slightly smaller than the 
entire boundary. 


6.2 Pair of pants 

For the pair of pants, we again study hnite temperature corrections by considering the departure 
of the mutual information between a region in one boundary, say boundary 3, and the whole of 
another boundary, say boundary 1, 


-^(^3 • 1) — Sa3 -f- iFi — Syigui, 


( 6 , 2 ) 


from the approximation suggested by (5.2). As for the two boundary case, we expect the 
main departure to come near transitions between different phases, where different geodesics are 
exchanging dominance in the calculation of the holographic entanglement entropies. For this case, 
the phase transitions depend on two parameters: the size of A^ and its location on boundary 
3. The different possible phases for S{A 3 ) and S{A 3 U Bi) are illustrated in hgures 10 and 11 
respectively. 

The calculation of the associated geodesic lengths can be easily automated using the description 
of the geodesic lengths as traces of corresponding SL{2,M.) group elements exploited in [27]. The 
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Figure 11: The geodesics giving the four possible phases of entanglement entropy of the union 
of a single interval Ai and the whole of boundary 2, in red, along with the event horizons added 
to satisfy the homology constraint, marked by dashed lines. 


lengths of the relevant geodesics can be found by computing the appropriate matrix products 
and traces. While the exact form of the answer is complicated and unilluminating, the general 
structure is fairly simple, being built mostly from polynomials in parameters encoding horizon 
lengths and the position of the interval. With a list of all contributing monomials in hand, hnding 
the length in the large L limit is equivalent to hnding the maximum of a set of linear functions. 
This calculation is implemented in Mathematica by performing a truncated series expansion. 

Of course, the full series can also be computed numerically. The results are summarized 
in hgure i2, which shows numerical results at hnite-temperature for deviations from the high- 
temperature approximation (5.2). The errors are indeed largest near the regions where nearby 
horizons are not exponentially close (i.e., where the causal shadows become large) and for intervals 
of size comparable to the thermal scale. Such regimes are close to phase transitions in the mutual 
information, where pairs of minimal curves exchange dominance. 

In addition to the bipartite entanglement, at large but hnite temperature one expects to hnd 
tripartite entanglement associated with the shadow region between the horizons. But as noted 
above the area of the pair-of-pants causal shadow is Acs = 27r in AdS units for all values of the 
moduli La- Chopping off the exponentially thin “arms,” it can be useful to model this region as 
an AdS-scale disk. This is quite reminiscent of the picture obtained in the tensor network model 
of the AdS vacuum in [23], where different spatial regions mostly had bipartite entanglement, 
with a residual multipartite component corresponding to an AdS-scale region. 
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Figure 12: Deviation of the mutual information I (A : Bi) between a region A = {03 G 
[0 — A(l)/2, 0 + A0/2]} in and the whole of boundary 1 from the piecewise-linear form implied 
by hgure 2 (lower right) for Li = 9, L 2 = 15, L 3 = 30. We plot the ratio between the error and 
the maximal mutual information (twice the entropy of Bi). Here 03 = 0 is the leftmost point in 
hgure 2 (lower right) and for comparison the inverse temperature 03 is ~ 1.3. The error is 
exponentially small in L, except in a region of thermal scale around certain phase transitions, 
where the order terms in (5.2) contribute. The diagonal lines with largest error occur where 
an endpoint of A leaves the region of B^ entangled with Bi. The vertical lines with similarly 
large error are along a plateau phase transition, as occurs in the two boundary case. 

7 Discussion 

Our main result is that, in the limit of large black hole horizons, the pair-of-pants wormhole 
in 2+1 gravity is dual to a CFT state formed by sewing together thermoheld doubles in one of 
the manners shown on the right of hgure 2 , or to the degenerate case that interpolates between 
them. We showed this by directly analyzing the CFT path integral dehning the state |S), and 
used consistency with bulk holographic calculations of the mutual information to argue that the 
E-wormhole dominates the associated bulk Euclidean path integral. We focused on the pair of 
pants for simplicity but - as will be discussed further below ~ it is easy to extend the central 
aspects of our discussion to more complicated wormhole spacetimes. 

We also focused on the case of circular boundaries, but the same conclusions apply to the 
planar case. In 2+1 bulk dimensions, such high-temperature n-boundary planar cases are just 
AdS 3 in non-standard coordinates corresponding to performing certain conformal transformations 
on the dual CFT vacuum that are singular at n points on the circle, with each segment running 
between two such singular points representing a distinct planar boundary. One may also consider 
wormholes having both planar and circular boundaries. 

Let us now briefly describe the extension to more general Riemann surfaces. Recall that a 
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general orientable Riemann surface S (other than the sphere or annulus) can be decomposed into 
pairs-of-pants. Let us think of E as the t = 0 slice of a wormhole spacetime with n boundaries 
each asymptotic to AdSs. Then the surface contains geodesics Hi (i = 1, 2,..., n) that define 
bifurcation surfaces of the event horizons for each boundary. In addition, it contains a number 
of internal geodesics. Each pair-of-pants decomposition of E corresponds to cutting E into 
pair-of-pants pieces along some set of these internal geodesics. It will be convenient for us to 
also cut along the Hi so that we in fact decompose E into n cylinders Q and some number of 
pair-of-pants pieces E/. In a somewhat redundant notation, we will refer to the geodesics forming 
the three boundaries of E/ as Hja for a = 1, 2, 3. Note that the set of Hia includes the horizons 
Hi. In this decomposition, the moduli space of the Riemann surface is parametrised by the 
lengths Lja of the Hja and the twists 9ia specifying the relative rotation between the two pairs of 
pants on the internal geodesics where we are sewing pairs of pants together. 

Each E/ has the same geometry as the causal shadow region lying between the three horizons 
in figure 6 as defined by the corresponding Lja- So each E± has area precisely 27r, independent of 
moduli. Any Hia which is long will lie exponentially close to another Hja (or another region of 
the same Hja) across the causal shadow region. As a result, a large number of such E/ can be 
sewn together without introducing an appreciable causal shadow or an appreciable reduction of 
the local energy density along each boundary. Away from the special points in each boundary 
corresponding to vertices in our previous discussion, the effective \TFD) temperature remains 
uniform in the round conformal frame specihed by the cylinders C*. Note that this is needed for 
consistency with the fact that the solution is precisely BTZ outside each horizon Hi, so that each 
boundary has constant energy density in our round frame^°. Some simple examples are shown in 
figures 13 and 14, the former being a 1-boundary wormhole whose causal shadow at t = 0 has the 
topology of a punctured torus. 

Interesting new behaviour can arise as a function of the twists as we sew together pairs of 
pants with the structure of the lower panel in figure 2, as we illustrate by example in figure 14. 
Consider for example a four-boundary wormhole with external horizons Hi, H 2 , H[, and split 
it into two pairs of pants along an internal geodesic H 3 = H'^ separating Hi ^2 from H[ 3 . We take 
L 3 ^ Li + L 2 , T 3 Hi + L 2 . If iLs and H'^ are identified via some twist 9, a given region in 
say Hi is entangled with a region in H^, which is in turn identified with some region in H'^. For 
generic 9 at large L 3 this will be entangled with some other region in iLg, which is then identified 
back to H^. For large L 3 we will cycle through the identification between H^ and H'^ many times 
before finally identifying the region with one of the other horizons {H 2 ,Hii,H 2 '). In the limit 
where L 3 is much larger than the external horizons, the identifications resulting from a general 
twist are chaotic and appear to give a fractal entanglement structure. It would be interesting to 
characterize these structures in more detail, and to relate this behavior to the well-known chaotic 

^°When the number of such pieces becomes comparable to the lengths Li of the horizons Hi, the qualitative 
effect on the CFT state jS) will depend on how these pieces of causal shadow are distributed along each boundary, 
and in particular on whether any regions of the boundaries do in fact remain far enough away to retain their 
\TFD) description. 
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Figure 13: A decomposition of the torus wormhole (left) with one boundary B into a single pair 
of pants El and a single cylinder Ci. Taking boundary 3 of Si to adjoin Ci, we see that there are 
three distinct moduli: the length L 3 of H, Li = L 2 = L, and a possible twist 6 . Tensor networks 
for 6 = 0 dual CFT states with large L3, L are also shown for L3 < 2L (middle) and L3 > 2L 
(right). In both cases, corresponding cyan and orange links are to be identihed as dictated by the 
twist angle 6 . For 6 = 0, this identihcation is reflection about the vertical axis through the center 
of each diagram. (Without this reflection, the spacetime is a punctured Klein bottle instead of a 
torus.) The cyan and orange links should be viewed as exponentially short, while the black links 
have length f3/2 set by the inverse effective temperature (3 of the black hole. So for 6 = 0 these 
identifications generate exponentially short closed loops which can be removed from the tensor 
without changing the state at leading order in large L3, L and central charge c. See discussion in 
main text below. 


dynamics of geodesics on compact hyperbolic spaces. 

Another snbtlety arises in cases like that shown in hgnre 14 a), where we consider the 
four-boundary system with zero twist, and take for simplicity Li = L[ and L 2 = L '2 with L 3 again 
very large. Then sewing together the two copies of hgnre 2 (lower right) indicates that Bi is 
entangled only with By and that B 2 is entangled only with B 2 '. As a result, taking A = BiU By 
and B = i ?2 Ui? 2 ', the CFT state has I{A : B) = 0 {ai leading order in large c and L). This result 
may seem is surprising from the bnlk point of view, as E contains a closed geodesic that rnns 
vertically around hgnre 14 (a), separating A from B. So HRT predicts I{A : B) = L/2G, where 
L is the length of this geodesic. This would be consistent with the above prediction as large L 3 
makes this geodesic exponentially short so that its length can be ignored at leading order. Note 
that this geodesic is short only for zero twist: we saw that for small-but-nonzero twist 6 a part 
of Bi is instead entangled in a local \TFD) state with part of B 2 ', so the mutual information is 
non-zero and grows as we scale up the La- Thus A and B can no longer be separated by a closed 
bnlk geodesic of negligible length^^. 


^^This is also clear from the fact that the local bulk geometry of these regions is essentially that of a segment of 
BTZ, and any such separating geodesic must traverse part of this segment of BTZ and thus have non-negligible 
length 
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Figure 14: (Left diagrams) Two pairs-of-pants S and E' are each cut along H^ y (red) and 
the pair of interior portions are sewn together along the cuts with twists 6 = 0,7 r/2, vr to form 
a four-boundary wormhole. Here we consider the case with > Li + L 2 , Ly > Ly Ly. In 
(a,c left) each pair of pants is bisected by an additional closed geodesic (not shown) that runs 
vertically around the diagram. In the high-temperature limit, the corresponding entanglement 
structure is given by identihying the outer boundaries H^ y of a pair of “eyeglass” diagrams, 
shown in red on the right-hand figure. This identification entails a twist 9 which is represented 
by the dot in each cut which are identified across the join in accordance with the twist. For 9 = 0 
(top) we infer that boundaries 1 and 2 are each entangled only with 1' and 2' respectively. For 
6* = I (middle) and the chosen values of La intervals within any given boundary are entangled 
with intervals in each of the others. The pattern of such entanglements become chaotic at generic 
9, though a twist in this setting never entangles two distinct intervals within the same boundary. 
For 9 = 71 (bottom) boundaries 1 and 2 are each entangled only with 2' and 1' respectively. 
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A - B --- C D — E. F . 

Figure 15: The entanglement strnctnre of the 4-bonndary wormhole (hgnre 14) with large 
L 3 = L 3 / for Li = Ly = L 2 = L 21 = \L^ and 9 = ^. The state is well-described by a tensor 
network analogons to hgnre 2 right. The state on any pair of bonndary intervals formatted 
in the same way (color, dots/dashes/solid lines) is a local piece of \TFD); the labels 1, l',2,2' 
indicate the boundaries connected to each TFD segment. The TFD intervals join at four vertices 
A, B, D, E located as shown. Each vertex connects the 3 local TFD states listed in the key below 
the diagram. C*, F are not vertices, but are simply in the middle of the indicated (relatively 
long) TFD intervals. Some TFD strips connect oppositely-oriented intervals, while some preserve 
orientation. 



Despite the above consistency, the appearance of such a short geodesic also suggests that 
our E-wormhole may not in fact be the dominant bulk saddle for the CFT state |E). It seems 
natural to conjecture that - unless forbidden by global features like a choice of spin structure - 
when E contains such exponentially short geodesics there will be another bulk saddle of smaller 
action where the geometry is modihed so that these geodesics are contractible when viewed as 
living on the boundary of the new saddle. That is, we conjecture that |E) in such cases is in fact 
dual to a bulk geometry with t = 0 surface E' built by cutting E along all exponentially small 
geodesics and capping off the resulting holes with small disks. This E'— “wormhole” may not then 
connect all the boundaries. From the tensor network point of view, the point is that the network 
obtained by gluing together two copies of hgnre 2 with no twist breaks up into two disconnected 
components, one connecting Bi and By and one connecting B 2 and i? 2 '- The remaining chains 
merely form closed loops. At leading order in large central charge c the properties of the state |E) 
are unchanged if we remove all chains that form closed loops rather than ending on boundaries. 
The resulting tensor network dehnes the manifold E'. The difference between |E) and |E') is 
then exponentially small at large c, and we conjecture the E'-“wormhole” to be the leading bulk 
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Figure 16: Another way to sew three \TFD) states together. Here the outer \TFD) has been 
sewn to itself at the ends as well as in the middle. The sewing at the ends creates features we call 
‘buds.’ Such buds are removed if one transforms the result to the round conformal frame. One 
may construct similar buds from the vacuum by applying a smooth conformal transformation 
approximating over some region the singular one that gives the infinitely long planar thermofield 
double state. 


saddle describing both states. This feature also arises for the punctured torus shown in figure 13 
for L3 < 2L (middle hgure) with vanishing twist 9. 

It is worth elaborating further on this last point. As noted in the caption for hgure 2, the 
diagrams in this paper include only a simple cartoon of the \TFD) tensor networks from e.g. [4]. 
The full tensor network for |E) obtained by sewing together \TFD) pieces as we describe will 
be correspondingly more complicated as well. In particular, returning to the simple example of 
two pairs of pants with very large L3 sewn together along the corresponding boundary, this full 
tensor network will certainly not factorize into unentangled states on BiBii and B 2 B 2 '. Instead, 
it will merely imply that the mutual information between BiBy and B 2 B 2 ' remains of order 1 at 
large central charge c. This is analogous to \TFD) below the Hawking-page transition where it 
describes two entangled copies of a thermal gas on pure global AdS3 backgrounds. Our conjecture 
is thus that the dominant bulk geometry at t = 0 is correctly predicted by removing parts of the 
full tensor network that fail to transmit mutual information of order c. We note that evaluating 
this criterion requires understanding the tensor structure of each node in the tensor network 
implied by the CFT dynamics; it is not apparent from the graph representation of the tensor 
network alone. 

So far we have considered tensor networks constructed by sewing together pair of pants 
networks in the way suggested by bulk wormhole geometries. But it is possible to consider a more 
general class of states dehned by sewing together high-temperature \TFD) states in arbitrary 
fashions. For example, one may sew a \TFD) to itself (or others) so as to introduce a ‘bud’ 
on the tensor network as shown in hgure 16. Second, some pieces of some \TFDYs - or even 
entire such states - may now be entirely internal to the tensor network, lengthening some chains 
and thus lowering the local temperature. In general, the chain length can then be non-uniform 
across any boundary. Together, these two effects recover the freedom to make arbitrary conformal 
transformations relative to the round conformal frame used above. That is, these more general 
states must be related to the states considered above rewritten in a more general conformal frame. 
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Finally, one may also generate non-orientable S by performing antipodal identifications on 
some circle boundary. For example, doing so one one boundary of a cylinder shows that the CFT 
state dual to the high-temperature AdSs geon (see e.g. [31]) is given by the thermofield double 
tensor network on a Mobius strip. The Mobius strip can of course be constructed by cutting 
open the cylinder along 0 = 0 and gluing the two ends back together with a half twist. It is an 
interesting question to what extent such gluing operations reproduce desired states when applied 
to particular e.g. MERA-like tensor network representations of states at finite temperature and 
finite central charge c. Results related to this issue will appear in [33]. 

While we have stressed the limit where all La become large, the discussion may be generalized 
to allow some La to remain small. The pair-of-pants CFT states |S) are then described by hgure 
2 (right) with the small-L boundaries contracted to points that merge with the vertices where 
the approximation by local TFDs breaks down. But regions of any large-L boundaries far from 
the new vertices remain well-described by the indicated local TFDs. One should be aware that, 
due to the possibility of bulk phase transitions like those described above, having some La small 
may make it less clear which bulk spacetime is in fact dual to |E). Nevertheless, the local TFD 
description of |S) remains valid. In particular, any entanglement of large-L boundaries with those 
having small-L will be conhned to intervals no longer than the effective thermal scale. The tensor 
network issue dual to uncertainties regarding bulk connectivity is that some new vertices may 
now be trivial in the sense that they no longer lead to order-c mutual information with the small 
boundary. When this occurs and creates a ‘bud’ as in hgure 16, the bud may again be absorbed 
into a neighboring vertex without changing the large-c structure of the state other than by acting 
with a conformal transformation^^. Similar comments apply to E having more boundaries or 
more general topology when some of the Lja remain small. 

Although we have discussed 2-|-1-dimensional bulk geometries above, but many of our 
considerations clearly apply to the higher dimensional case as well. In particular, sewing together 
high-temperature |TL’D)’s dehnes a zoo of interesting states |E) and conformal geometries E. 
And it is again natural to conjecture the CFT states |E) dehned by such sewing operations to 
be dual to E'-wormholes dehned by having a moment of time-rehection symmetry on which the 
induced geometry differs from (planar) Schwarzschild-AdS only by small corrections outside a 
hnite number of AdS-scale regions. But much remains to be understood and the details will prove 
interesting to explore. In particular, one would like to hnd an algorithm that takes the tensor 
network naturally associated with |E) dehned by the above gluing procedure and turns it into 
one in which the geometry of E' is manifest - e.g., with the tensor network providing a cellular 
decomposition of E' in terms of AdS-scale cells [34-37]. One wonders if solving the Euclidean 
Einstein equations to construct E' from E can be related to a renormalization-group how on 


^^This creates a local connection between antipodal points on the boundary. The fact that the bulk geodesic 
between antipodal points is short in the large temperature limit can be seen from the explicit formula for the 
geodesic lengths in [32], although it is incorrectly stated there that the length of the geodesic through the 
identification is always longer than the one outside the horizon. 

^^It is then the diagram without the bud that describes the round conformal frame. 
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tensor networks akin to those discussed in [38, 39]. 
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A The horizons Hi^H2 in BTZ coordinates 

We now compute the parameters characterising the region S+, corresponding to half of the f = 0 
slice of the three-boundary wormhole with horizons Hi, H 2 , H 3 , of respective lengths Li, L 2 , T 3 . 
In the main text we ordered the lengths so that T 3 > Li, L 2 ; this assumption is relaxed here. 
The region is bounded by three geodesics Gab, running between the boundary components 
labelled by a and b, and meeting horizons Ha, Hf, orthogonally. 

We use the metric (3.1) with H 3 lying at p = 0 and (^ 13,^23 lying at x = ±^. Thus x G 
Consider a geodesic parameterised by arclength s, using a dot to denote differentiation 
with respect to s. From translation invariance, there is a conserved quantity (1 -|- p^)x, which for 
geodesics with both endpoints at p = 00 is given by -y/I + pg, where po > 0 is the minimal value 
of p. The geodesic is then given by 


P 


Po coshs. 


X = Xq + tanh ^ 


tanh s \ 
+ Po) 


(A.l) 


Consider hrst the geodesics G 12 and that corresponding to (see hg. 17 for the various 
relevant geodesics and quantities). The endpoints of G 12 lie at x = Xi, X 2 , and it will be convenient 
to parametrise these by the centre x = and the half-width Ax = . We intend to hnd 

Li in terms of these parameters, and along the way will also obtain the minimal distance di 3 
between and H^^, as well as the position at which G 12 and i7i+ intersect. 

The geodesics are given by 


Gi 2 ■ p — 


cosh s 
sinh Ax ’ 


X = X -f tanh ^ (tanh Ax tanh s ), 


s e 


H 


1 + : p = sinh dig cosh; 


X = 


, _i f tanhs 
+ tanh —^—— 

y cosn di 3 ^ 


Li 

0<.<y, 


(A.2) 

(A.3) 
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Figure 17: The geodesics Gab bounding the patch in BTZ coordinates. The horizons Ha+ 
are also shown. The p direction runs vertically, and x horizontally. The positions where the 
geodesics intersect are labelled with the arclength along each curve, measured from the deepest 
point (minimal p), and di 3 ,d 23 mark the minimal distances between the horizons. 


with the constraint that they intersect at right angles at the endpoint of i7i+, where the arclength 
along Hi+ is s = and along G 12 is s = si, say, where si < 0 . 

The condition that two geodesics intersect orthogonally determines the value of p at which 
they meet in terms of the conserved quantities for the two geodesics; for Hi^ and G 12 h gives 

p^ = sinh^ di 3 + coth^ Ax = cosh^ ^13 + csch^ Ax at intersection. (A.4) 


We now get two equations from identifying the value of p at intersection with the values of p for 
Gi 2 at s = Si, and for Hi^ at s = ^. A third comes from identifying the x coordinates at these 
same arclengths: 


sinhsi = 


nhdiq sinh — 

id 2 

■sinh Ax cosh(ii 3 . 


X + tanh ^(tanhAxtanhsi) = —— + tanh 


-1 



(A.5) 

<^13) 

(A. 6 ) 

^ tanh ^ ^ 


1 cosh (ii 3 j 

(A.7) 


We then solve for Li, di 3 , si in terms of x. Ax to obtain 


Li sinh (¥ + x) cosh Ax 

cosh— = - ^tanhai 3 =- yj -tanhsi = 

2 sinh Ax ’ cosh(^+x) 


tanh Ax 


tanh (^ + x) 


(A.8) 


It is straightforward to translate these results into expressions for L 2 , the distance ^23 between 
horizons if 2 + and ii 3 +, and S 2 , the arclength along G 12 at which it intersects ii 2 +: 


cosh 


L2 


sinh 


X I 

. V. -, tanhd23 

sinh Ax 


cosh Ax 
cosh — xY 


tanh 


S2 


tanh Ax 
tanh — x) 


(A.9) 
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Finally, the above can be inverted to find x, Ax, and dab (where, in particular, du = S 2 — si 
is the minimal distance between Hi^,H 2 +) in terms of Li,L 2 ,L 3 . 


sinh X = 


sinh Ax = 


cosh(ii2 = 


(cosh ^ — cosh sinh ^ 
3 sh^ ^ + cosh^ ^ 


^ + 2 cosh ^ cosh ^ cosh ^ 
sinh ^ 


+ cosh^ 


r - 


(A.IO) 

(AAl) 

(AA2) 


In particular, the explicit inversion shows that the mapping between (xi,X 2 ) and (Li,L 2 ) is 
bijective and smooth. 

We may now work out the asymptotic values of these quantities in the limit where all lengths 
La are large. The typical expressions reduce to sums of exponentials of linear combinations of La, 
so there are separate regimes depending on the relative sizes of the exponents; these turn out 
to be three regimes where one horizon is longer than the sum of the others (Li > L 2 + L^ and 
permutations), and the regime where no horizon is dominant in this way. 


2 ^2 Li> L2 + T3 

-^ + |exp(-^ 

L2>Li + L 3 

L 1 —L 2 

4 

otherwise 

exp 

Li > L2 + L3 

exp( 

L2 > Li L3 

L,-Li-L2 +log2 

L3 > Li + L2 

exp( 

otherwise 

L 3 -L 1 -L 2 +2 log 2 

L3 > Li -\- L2 

2exp(-y) 

Li > L2 -\- L3 

2exp(-y) 

L2 > Ti + L3 

2 exp( 

otherwise 


(A.13) 


(A.14) 


(A.15) 


The corrections in each case are exponentially small in the La, except when L 3 — Li — L 2 is order 
one, for example. 

The interval [xi,X 2 ] looks qualitatively different in each of the four regimes. When L 3 > 
Li + L 2 , it is long (the same order as the horizon lengths), and at a generic position. When 
Li > L 2 + T 3 , it is exponentially short, and also close to the right end of the strip; it is similarly 
short and close to the left end when L 2 > L 1 +L 3 . In the remaining regime, it is again exponentially 
short, but in a generic position. 
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